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1.  Introduction 

In  a  previous  paper  IlO]  we  studied  the  problem  of  determining  the  distribution 
of  the  lighted  portion  (measure  of  vacancy)  of  the  circumference  of  a  circle 
having  a  source  of  light  at  the  center  and  shadowing  arcs  generated  by  randomly 
placed  disks  within  the  circle.  The  centers  of  the  disks  are  assumed  to  be 
uniformly  distributed  on  the  plane  and  their  diameters  are  i.l.d.  random  vari¬ 
ables  having  a  general  distribution  independent  of  the  location.  Formulae  were 
given  for  the  numerical  determination  of  all  moments  of  the  measure  of  vacancy. 

The  present  paper  generalizes  the  methodology  of  obtaining  moments  of  the  measure 
of  vacancy  to  cases  where  a  more  general  type  of  random  objects  are  casting 
shadows  on  star-shapped  curves  on  the  plane.  He  start  the  development  in  the 
next  section  by  defining  a  class  of  random  objects  in  the  plane,  which  satisfy 
certain  initial  conditions.  In  particular,  we  consider  the  class  of  Poisson 
random  objects  for  which  the  methodology  of  the  present  paper  is  developed.  He 
introduce  also  in  this  section  the  basic  concepts  related  to  the  shadowing  pro¬ 
cess.  In  section  3  we  derive  general  formulae  for  the  probability  of  the 
simultaneous  visibility  of  n  points  in  the  plane.  These  formulae  are  applied 
in  section  4  to  the  development  of  recursive  formula  for  the  moments  of  measures 
of  visibility  of  curves  in  the  plane.  Examples  of  the  application  of  the  method¬ 
ology  presented  in  the  report,  different  from  the  one  worked  out  In  [10]  will 
be  given  in  following  reports. 

The  literature  on  shadowing  processes  is  quite  limited.  Chernoff  and  Daly 
[2]  studied  the  distribution  of  length  of  shadows  of  disks  on  a  line.  The 
shadowing  problem,  however,  is  a  special  case  of  the  general  coverage  problem  on 
which  there  is  extensive  literature.  In  particular  we  refer  to  the  studies  of 
Robbins  [5,6],  Ailam  [1],  Siegel  [7,8]  and  the  monograph  of  Solomon  [9]  which 
summarizes  many  of  the  important  results. 


2.  Random  Objects  and  Their  Shadows 


Consider  a  countable  set  of  disks,  S,  scattered  on  the  plane.  Each  disk 
is  specified  by  the  polar  coordinates  (p,0),  0£p<®,  -n_<0<fl,  of  its  center  with 
respect  to  an  origin  0  (observation  point,  source  of  light);  and  by  its  diameter 
y,  0<y<a,  (a<°°) .  Let  Rq  denote  the  collection  of  Borel  subsets  of  the  sample 
space 

(2.1)  Rq  *  ((p.O.y):  0<p«»,  -*£0<n,  0<y<a}. 

For  every  BeRq  denote  by  N{B)  the  number  of  disks  such  that  (p,6,y)eB  (satisfy¬ 
ing  condition  B) .  If,  for  every  BeR^,  N{B}  is  a  random  variable,  then  the 
elements  of  S  are  called  random  disks.  In  particular,  the  disks  are  called 
Poisson  random  if  for  every  BeRq,  N { B }  has  a  Poisson  distribution  with  mean 

(2.2)  v(B)  -  u///G(dy|p,0)H(dp,d0) 

B 

where  p,  0<p«»>,  is  an  intensity  parameter;  H(p,0)  and  G(y|p,0)  are  the  sigma-finite 
measure  of  the  locati-n  coordinates  (p,0)  and  the  conditional  c.d.f.  of  the 
diameter,  y,  given  (p,0),  respectively. 

The  special  case  in  which  H(dp,d0)  -  pdpd0  and  G(y|p,0)  «  G(y),  is  called 
the  standard  case.  The  standard  case  is  the  one  in  which  the  centers  of  the 
disks  are  uniformly  distributed  on  the  plane  and  their  diameters  are  independent 
of  their  location.  In  this  case 

(2.3)  v(B)  -  u/8  HB(y)dG(y), 

where  w  is  the  mean  number  of  disks  per  unit  area;  K^(y)  is  the  area  of  the 
region  in  which  disks  of  diameter  y,  satisfying  condition  B,  are  centered. 

The  development  of  the  theoretical  framework  is  not  restricted  to  random 
disks  in  the  plane  but  can  be  generalized  to  Poisson  random  objects  in  the  plane, 
or  e'  n  more  generally  in  the  space,  by  replacing  the  parametric  vectors  (p,6ty) 
to  a  general  parametric  vector  (nj,...,n  ),  which  characterise  the  location. 
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orientation,  shape  and  size  of  a  family  of  objects.  In  such  cases  Rg  and  Rg 
are  properly  generalized.  In  order  to  simplify  the  presentation,  we  restrict 
attention  in  the  present  study  to  the  family  of  random  disks  in  the  plane.  Most 
of  the  formulae,  however,  are  valid  for  general  families. 

A  natural  requirement  for  shadowing  processes  is  that  the  source  of  light 
(the  origin)  is  uncovered.  We  therefore  introduce  the  initial  condition 

(2.4)  CQ  =  {(p,0,y);  |lP<c°»  »  0<y<a} 


and  assume  that  all  random  disks  satisfy  Cg.  In  special  cases  one  considers 
more  stringent  initial  conditions,  specified  by  sets  C  contained  in  Cg.  For 
example,  in  the  previous  paper  of  Yadin  and  Zacks  [10]  the  initial  condition 
considered  is 

(2.5)  c  *  {(p,e,y);  |<p<i-^;  -n<e<ir,  o<y£i)  . 


Let  P  be  a  point  in  the  plane.  P  is  said  to  be  visible  (in  light)  if  the 
line  segment  OP  does  not  interset  any  random  disk:  The  set  of  all  visible  points  in 
direction  s,  -it<s<ir,  starting  at  the  origin,  is  called  a  line  of  sight.  Lg.  Let 
P=(r,s)  be  a  point  in  the  plane,  in  direction  s  and  distance  r  from  the  origin. 

A  random  disk  (p,8,y)  intersects  the  line  segment  OP  if,  and  only  if,  (p,6,y) 
belongs  to 

(2.6)  B(r,s)  *  {(p,e,y);  0<p<p(6,y,r,s) , 

-iT^eor,  0<y<a}  , 


where 


r  cos(e-s)  +  ((£)  -r2sin2(6-s))'*,  if  |e-s|<_tan  1(^~) 


(2.7) 


p(0,y,r,s)  -  2  sin  | 0-s 


2r 

,  if  tan  1(2— )< |e-s|jt  j 


1 

2 
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Accordingly,  a  point  P*(r,s)  is  visible  if  N{B(r,s)nC}  =  0.  Hence,  a  line  of 
sight  Lg  has  magnitude 

(2.8)  ||L  ||  *  sup{r;  N{B(r,s)nC)  -  0}  . 

s 

Notice  that  j | L  ||  is  a  random  variable, 
s 

Consider  a  curve,  C,  in  the  plane  such  that  each  ray  originating  at  0 
intersects  C  at  most  at  one  point  (star-shaped  curve).  Such  a  curve  is  specified 
by  a  continuous  positive  function  r(s) ,  s^s^s^.,  -it<sL<Sjj<ir,  where  (r(s),s)  are 
the  polar  coordinate  of  points  on  C.  We  further  consider  smooth  curves,  that 
is  r(s)  is  almost  everywhere  differentiable.  We  assume  that  the  initial  con¬ 
dition  C  in  Cq  implies  that  C  is  uncovered  by  random  disks.  Define  the  indicator 
function  I(s),  s^^s^Sy  so  that  I(S)=1  if  (r(s),s)  is  a  visible  point  and  l(s)=0 
otherwise.  The  total  measure  of  the  visible  part  of  C  is  given  by 

(2.9)  V{C}  -  /U  I(s)(r2(s)  +  (r’(s))2)%ds 

SL 

V{C}  is  the  visibility  measure  of  C.  In  coverage  problems  it  is  known  also  as 
the  measure  of  vacancy  (Ailam  {1 ] ) . 

3.  Visibility  Probabilities  for  Poisson  Random  Objects 

Let  P=(r,s)  be  a  point  in  the  plane.  Under  the  Poisson  randomness  assump¬ 
tion,  the  probability  that  P  is  visible  is 

(3.1)  Q(r,s)  -  P[N{B(r ,s)nC}  -  0] 

-  exp{-v(B(r,s)nC))  . 

From  this  function  one  can  immediately  obtain  the  distribution  of  ||Lg||. 

Indeed , 

(3.2)  P{  |  |  Lg|  |  >1}  ■  Q(£,s)  ,  0 
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Furthermore,  according  to  (2.2)  and  (2.6)  one  obtains  that 

a  tt/2  p(0,y,r,s) 

(3.3)  v{B(r,s)nCn)  *  /  J  /  G(dy|p,0)H(dp,d0) . 

U  0  -t/2  y/2 

In  particular,  for  the  standard  case 

(3.4)  v{B(r,s)nC0>  *  p£r  , 

where  £;  is  the  expected  diameter  of  random  disks.  This  implies  that  in  the  standard 
case  the  distribution  of  ||(-sl|  is  exponential  with  mean  1/p,  where  p  ■  p£. 

This  special  result  is  well  known  and  may  be  obtained  directly  from  the  Markov¬ 
ian  nature  of  the  phenomenon  (see  Feller  [3,  p.  10]). 

Let  C  *  {(r(s),s),  s^fs^s^}  be  a  specified  smooth  curve  in  the  plane. 

Let  s^  and  s^,  s^^s^s^s^,  be  the  orientation  coordinates  of  two  specified 
points,  and  P^,  on  the  curve.  In  the  following  development  we  consider  only 
initial  conditions  C  which  specify  the  set  of  all  (p,0,y)  points  which  may  cast 
shadow  on  the  curve  C.  Note  that  v{C}«*>  for  all  curves  of  finite  length.  A 
necessary  and  sufficient  condition  for  P.^  and  to  be  simultaneously  visible 
is  that 

N{(B(r(s1),s1)uB(r(s2),s2))nC)  ■  0  . 

Since  B(r(s^),s^)  and  B(r(s2>,s2)  are  generally  not  disjoint  it  is  more  con¬ 
venient  to  consider  a  union  of  disjoint  sets  obtained  in  the  following  manner. 

Define,  for  an  s  in  (s  ,s..], 

L»  U 

(3.6)  B+(s)  -  {(p,e,y):  Js^feTg7  <p<®;  0<y<a} 

and 

(3.7)  B_(s)  -  { (p,0,y) :  2sin(s-0)'  <p<,0;  8L<6<8*  0<y<«}  . 
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B+(s)  and  B  (s)  are  the  sets  of  all  possible  random  disks  on  the  right  and  on 
the  left  of  a  ray  with  orientation  s,  which  do  not  intersect  it.f  The  set  under  con 
sideration  is  the  complement  with  respect  to  C  of  the  union  of  the  disjoint  setB 
B  (s^)uB+(s2),  B+(s^)nB_(s2) .  Accordingly,  the  probability  that  and  are 
simultaneously  visible  is 
(3.8) 

P(s1,s2>  *  exp{-[v{C)  -  v{(B_(s1)uB+(s2))nC)  -  v{(B+(s1)nB_(s2))nC}]}  . 

Generally,  if  sL<s^< . . . <sn<s^  are  the  orientation  coordinates  of  n  points  on  the 
curve  C,  the  probability  that  they  are  simultaneously  visible  is 

P(s.,...,s  )  =  exp{-[v{C)  -  v{ (B  (s.)uB  (s  ))nC) 
in  -  1  +  n 

n-1 

-  E  v{(B+(si)nB_(s1+1))nC)]}. 

4.  Moments  of  Visibility  Measures 

According  to  (2.9)  the  expected  values  of  the  visibility  measure,  V{C),  Is 

s, 

(4.1)  E{  V{C} }  «/U  E{I(s)  K(s)ds, 

8L 

2  2  k 

where  l(s)ds  ■  [r  (s)  +  (r'(s))  ] *ds  is  the  inflnitesmal  length  of  C  at  s  and 
E(I(s)}*  Q(r(s),s).  It  may  often  be  more  convenient  to  compute  Q(r(s),s) 
according  to  (3.8),  in  which  we  substitute  s^*s2-s.  In  this  case, 

(4.2)  Q(r(s),s)  -  P(s,s) 

■  exp{-[v{C)  -  v{(B_(s)uB+(s))nC}]). 

Generally,  for  every  n_>2  the  n-th  moment  of  V{C}  is  given  by 


(4.3) 


n!  / 


{8L<S1<...<8nlSu} 


/  P(s, ,...,s  )  n  £(s.)ds. 

1  n  i-1  1  i 


Indeed, 

visible. 


n 

E{  n  I(s  )}  is  the  probability  that  all  the  n  points  are  simultaneously 
i-1  1 

According  to  (3.9), 


(4.4) 


-v{C}_ 


n-1 


P(3l,...,sn)  -  e  ^i’W’ 


where 


^(sl*sn)  =  exP(vHB_(s1)uB+(sn))nC}} 

(4.5)  and 

i)j(s,t)  =  exp{v{(B+(s)nB_(t))nC}} 


Let  Q{C}  denote  the  probability  that  the  whole  curve  C  is  visible.  Notice  that 
Q{C}  -  exp{-v{C}} .  From  (4.3)  -  (4.5)  we  obtain 


s  s 

E(Vn{C}}  -  n!Q{C)/U  dsn  JU  dsx  •  *(8^ ^(s^)/1*  ds2 

SL  SL  81 

<4’6> 

•  i(s2)i(i(a1,s2)  /n  .  .  .  Jn  dsn_2 

s2  sn-3 

*  £<sn-2),J’(sn-3’Sn-2)/n  d8n-l£(8n-l)*(sn-2*Vl)£(sn)+(Vl>Sn) 

sn-2 

Define  recursively  the  functions 

n1(s,t)  -  £(s)if»(s,t)  »  SjlsiCiSu 

(4.7) 

^(s.t)  -  Jtn1(s,v)nk_1(v,t)dv. 
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Combining  (4.6)  and  (4.7)  we  obtain,  for  every  n^2, 

s  s 

(4.8)  EtV^CH  =  n!Q{C}JU  JU£(sK(s,t)n  1(s,t)ds  dt 

sL  S 

In  cases  whter  B  (s)nC  and  B+(t)nC  are  disjoint  sets  for  every  s<t  one  can 
derive  simpler  expressions  for  the  moments. 

Finally,  let  1{C }  denote  the  length  of  C.  The  distribution  of  V{C}/£{C) 
is  concentrated  on  [0,1],  with  jumps  at  the  two  end  points  0  and  1  and  absolutely 
continuous  elsewhere.  It  follows  that 

(4.9)  lim  E {V^CH/AC}  =  P{V{C}  =  £{C}}  =  Q{C) 
n-*» 
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